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ABSTRACT 

We use hydrodynamical simulations of different dark energy cosmologies to investigate the 
concentration-mass (c — M) relation in galaxy clusters. In particular, we consider a reference 
ACDM model, two quintessence models with inverse power-law potentials (RP and SUGRA), 
and two extended quintessence models, one with positive and one with negative coupling with 
gravity (EQp and EQn respectively). All the models are normalized in order to match CMB 
data from WMAP3. We fit both the dark matter only and the total mass profile with a NFW 
profile, and recover the concentration of each halo from the fit using different definition. We 
consider both the complete catalog of clusters and groups and subsamples of objects at dif- 
ferent level of relaxation. We find that the definition itself of the concentration can lead to 
differences up to 20% in its value and that these differences are smaller when more relaxed 
objects are considered. The c—M relation of our reference ACDM model is in good agreement 
with the results in literature, and relaxed objects have a higher normalization and a shallower 
slope with respect to the complete sample. The inclusion of baryon physics is found to influ- 
ence more high-mass systems than low-mass ones, due to a higher concentration of baryons 
in the inner regions of massive halos. For the different dark energy models, we find that for 
ACDM, RP and SUGRA the normalization of the c — M relation is linked to the growth 
factor, with models having a higher value of <7sD + having also a higher normalization. This 
simple scheme is no longer valid for EQp and EQn because in these models a time dependent 
effective gravitational interaction, whose redshift evolution depends on the sign of the cou- 
pling, is present. This leads to a decrease (increase) of the expected normalization in the EQp 
(EQn) model. This result shows a direct manifestation of the coupling between gravity and 
the quintessence scalar field characterizing EQ models that can be in principle investigated 
through the analysis of the c — M relation. 
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1 INTRODUCTION 

The internal properties of dark matter halos are known to 
reflect their formation history and thus the evoluti on of the 
background cosmology. iNavarro. Frenk & Whitel d 1996k (hereafter 
NFW) found that the dark matter profile of a halo can be character- 
ized by a scale radius, which is linked to the virial radius through 
the concentration of the object. The concentration of a dark matter 
halo is related to the mean density of the universe at the halo for- 
mation time. Moreover, the presence of baryons influences these 
internal properties because they tend to concentrate in the inner re- 
gions, increasing thus the concentration of the halo. 

Because of the hierarchical nature of structure formation and 
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the fact that collapsed objects retain information on t he background 
avera ge matter density at the time of their formation (INavarro et al.l 
1996), concentration and mass of a dark matter halo are re- 
lated. Since low-mass objects form earlier than high-mass ones, 
and since in the past the background average matter density was 
higher, low-mass halos are expected to have a higher concentra- 
tion compared to high-mass ones. For a given mass, the evolution 
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the results of iV-body numerical simulations which find, at z = 0, 
a concentration-mass relation c(M) oc M°, with a ~ —0.1 
dDolag et alj|2004l ; iGao etai]|2008l ; Izhao et alJl2010T> . with a log- 
normal scatter ranging from 0.15 for relaxed systems to 0.30 for 
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disturbed ones Jjin 3 l2000h . This is true for idealized, dark mat- 
ter only halos, but in real objects the impact of baryon physics is 
strong, in particular in the inner regions. Indeed, even if the ICM 
is described by a cored /3-model profile, in the central region it 
tends to cool and to form stars, which accumulate in the core of 
the cluster, thus increasing the total concentration. Moreover, it 
is supposed that this baryon contraction does also affect the dark 
matter component, through an adiabatic contraction effect, leading 
thus to an increase of the dark matter con centration. The adiabatic 
contraction model has been described in iBlumenthalet all d 19861) 
andlRvden & GunrJ ( fl987h . but, for an updated treatment, see also 
lFedelild2012f> ! 

TV-body simulations have been carried out by several authors 
in order to study the c — M rel ation in dark matter halos with sizes 
of galaxy groups and clusters. [Dolag et ail (2004) performed sim- 
ulations with different cosmological models in order to verify the 
effects of dark energy dynamics. For ACDM they found that the 
dependence of concentration on mass and redshift can be fitted by 
a relation of the form (1 + z)c(M) = M a , with a « —0.1. For 
dark energy models, they found that the halo concentration depends 
on the dark energy equation of state through the linear growth fac- 
tor at the cluster formation redshift, D+iz^nii). iNeto et al.1 d2007h 
also noted that non-relaxed objects have a lower concentration 
and a higher scatter with respect to relaxed ones. Macc io et al.1 
d2008h made a comparison between concentrations in the WMAP 1 , 
WMAP3 and WMAP5 cosmologies in order to study the effects of 
different cosmological para meters (in particular <j$) on the c — M 
relation. iDuffv et aU ((2008) found that the concentrations of their 
halos is lower than the one inferred from X-ray observations and 
addressed that fact to the effect of baryon physics that was miss- 
ing in their simu lations. Actually , despi te an initial compatibility 
of the results in IPratt & Arnaudl ( 2005h a nd Pointecouteau et al] 
(2005) with the predictions of iDolag et al.l d2004l) . considering the 
most recent data sets there is still poor agreement between the ob- 
served c — M relation and the predicted one, with the first hav- 
ing, in general, a steeper slope and a higher normalization com- 
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the impact of baryon physics on the structure of dark matter halos, 
and to try to reconcile the discrepancy between the expectations 
from dark matter only simulations and the concentrations i nferred 
from observations, were carried out in buffv et al l (|2010|) . How- 
ever, they showed that even including baryon physics in the sim- 
ulations they cannot reproduce both observed concentrations and 
stellar fraction in galaxy groups and clusters. In general, including 
baryons they found lower concentrations with respect to the dark 
matter only case. The triaxiality of the matter distribution and the 
presence of substructures within the host halo virial radius could 
bias the estimates of the mass and concentration obtained through 
lensing techniques at the level of few per cent, but cannot rec oncile 
the observed and simulated relations (see lGiocoli et al.ll2012l) . 

Since the concentration of a halo is linked to the background 
density of the universe at the time it collapsed, and since different 
dark energy models predict different evolutions of the cosmolog- 
ical background, it is interesting to investigate the impact of dark 
energy on the c — M relation. Moreover, since some dark energy 
models can also affect the linear and non linear evolution of the 
density fluctuations, leaving some imprints in collapsed structure, 
one can think about using the c — M relation as a cosmological 
probe, orthogonal to others that are commonly used. TV-body cos- 
mological simulations of extended quintessence models, including 



the effects on the c—M relation, were presented in lLi et ai1d201 ll) . 
iBasilakos etal] ( 2009 ) performed a theoretical study on the spher- 
ical collapse model beyond ACDM and found that the observed 
c—M relation could indicate the possibility for a clustered dark 
energy. 

In this paper we address this issue by studying the c — M 
relation in galaxy clusters extracted from dark matter only and hy- 
dr odynamical simulatio ns of five dark energy models introduced 
in lDe Boni et alj d201 ll) (hereafter Paper I). These simulations al- 
low us to evaluate at the same times the effects of baryon physics 
and dynamics of dark energy on the internal properties of galaxy 
clusters. 

The paper is organized as follows. Dark energy models and 
numerical simulations are presented in Sect. [2] After introducing 
different methods to define the concentration in Sect. [3] along with 
a non-parametric definition in Sect. [4] we compare our results for 
the ACDM model with other results from literature in Sect. [5] For 
the different dark energy cosmologies, we analyse the dark matter 
only runs in Sect. [6] and the hydrodynamical runs in Sect. [7] We 
discuss the results in Sect.[8]and draw our conclusions in Sect. [9] 



2 THE COSMOLOGICAL MODELS 

We consider the same cosmological models discussed in Paper I. 
Here we recall only the main features of the different models, and 
refer to Paper I for more details. 

As a reference model we us e the concordance A CDM model, 
adapted to the WMAP3 values dSpergel et ai]|2007l) , with the fol- 
lowing cosmological parameters: 

• matter density: Qom = 0.268 

• dark energy density: £Iqa = 0.732 

• baryon density: Oot = 0.044 

• Hubble parameter: h — 0.704 

• power spectrum normalization: erg = 0.776 

• spectral index: n 3 = 0.947 

This model is characterized by the presence of a dark energy 
component given by a cosmological constant A, with a constant 

WA = — 1. 

The second case is a model with dynamical dark energy, given 
by a quintessence scalar field <j> with an equation of state w = w(a) 
dWetterich 1988: iRatra & Peebleslll988l) . As in Paper I, as poten- 
tials for minimally coupled quintessence models, we consider an 
inverse power-law potential 
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the so called RP potential dRatra & Peebles! 1 19881) . as well 
as its generalizatio n suggested by supergravity arguments 
dBrax & MartirJll999h . known as SUGRA potential, given by 



V(4>) = 



M 



■ exp(47rG0 
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where in both cases M and a ^ are free parameters (see Table [T] 
for details). 

The third possibility we cons ider is the case in which <j> inter- 
acts n on minimally with gravity (Wetterich 1988; Boisse au et al.1 
2000). In particular we refer to the extended quintessence (EQ) 
models described in IPerrotta et all d2000h . IPettorino etail d2005t) 
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Figure 1. Correction to the gravity constant for the two extended 
quintessence models, EQp (cyan) and EQn (red), as expressed in equation 



and IPettorino & B accigalupil d2008h . The parameter £ represents 
the "strength" of the coupling (see Table Q] for details). In partic- 
ular we consider here a model with positive coupling £ > (EQp) 
and one with negative £ < (E Qn). For an extensive line a r treat - 
ment of EQ models we refer to IPettorino & Baccigalupil (2008). 
Here we only recall for convenience that EQ models behave like 
minimally coupled quintessence theories in which, however, a time 
dependent effective gravitational interaction is present. In particu- 
lar, in the Newtonian limit, the gravitational parameter is redefined 
as 
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Here the coupling F(<j>) is chosen to be 



F(<t>) = - + ar 
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with k = %txG* , where G* represents th e "bare" gravitational con- 
stant (Espo sito-Farese & Polarski 200l|). 

For small values of the coupling, that is to say ( < 1, the latter 
expression becomes 
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which manifestly depends on the sign of the coupling £. We note 
that, since the derivative of the RP potential in equation Q} with 
respect to tj> is dV(4>)/d<f> < 0, we have <j) 2 < This leads to the 
behaviour of G/G* shown in Fig.Q] Note that the corrections are 
only within the percent level. 

For the ACDM model, the linear density contrast of a top-hat 
spherical collapse model, S c , converges to the Einstein-de Sitter 
value 1.686 at high redshift, but it is lower than the Einstein-de 
Sitter value at low redshift due to the presence of the cosmological 
constant which starts to dominate. For minimally coupled and 
extended quintessence models, the value of 8 C at low redshift is 
lower than the ACDM one, depending on the model. In any case, 
at z = these differences are very small, well below 1%. For an 



Table 1. Parameters for the different cosmological models: a is the expo- 
nent of the inverse power-law potential; £ is the coupling in the extended 
quintessence models; wq is the present value of the equation of state pa- 
rameter for dark energy; erg is the normalization of the power spectrum. 



Model 


a 




w 




ACDM 






-1.0 


0.776 


RP 


0.347 




-0.9 


0.746 


SUGRA 


2.259 




-0.9 


0.686 


EQp 


0.229 


+0.085 


-0.9 


0.748 


EQn 


0.435 


-0.072 


-0.9 


0.729 



extensive treatment of EQ models see Pace et al. (in prep.). 

We trimmed the parameters of the four dynamical dark energy 
models in such a way that Wo ~ w(0) ~ —0.9 is the highest value 
still consistent with present observational constraints in order to 
amplify the effects of dark energy. The parameters Slom, Qoa, fiob. 
h, and n s are the same for all models, but since we normalize the 
power spectrum to CMB data from WMAP3, this leads to different 
values of as for the different cosmologies: 
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assuming zcmb = 1089, where the growth factor D + is normal- 
ized to unity today. In principle, in EQ models there could be a 
scale-dependent linear growth rate which is neglected in equation 
((6). Given the values of the coupling, we do expect this dependence 
to be small. Table [TJ lists the parameters chosen for the different 
cosmological models. 

In order to study the formation and evolution of large 
scale structures in these different cosmological scenarios we use 
A-body and A-body + h ydrodynamical simul a tions performe d 
with the GADGET- 3 code ISpr ingel et ail 1200 ll ; ISpringell 120051) . 
which makes use of the en tropy-conserving formulation of SPH 
JSpringel & Hemquisdl2002l) . The hydrodynamical simulations in- 
clude radiative c ooling, heating by a un iform redshift-dependent 
UV background jHaardt & Madaull 19961) . and a treatment of star 
formation and feedback processes. For the dark matter only simu- 
lations, we simulated a cosmological box of size (300 Mpc /i -1 ) 3 , 
resolved with (768) 3 dark matter particles with a mass of rridm ~ 
4.4 x 10 9 Mq For the hydrodynamical simulations, we con- 
sidered a cosmological box of size (300 Mpc h' 1 ) 3 , resolved 
with (768) 3 dark matter particles with a mass of nidm ~ 3.7 x 
10 9 Mq and the same amount of gas particles, having a mass 
of m ga3 « 7.3 x 10 8 M h~\ 

Using the outputs of simulations, we extract galaxy clusters 
from the cosmological boxes, using the spherical overdensity cri- 
terion to define the collapsed structures. We take as halo cen- 
tre the position of the most bound particle. Around this particle, 
we construct spherical shells of matter and stop when the total 
(i.e. dark matter plus gas plus stars) overdensity drops below 200 
times the mean (as opposed to critical) background density de- 
fined by Q. m poc\ the radius so defined is denoted with i?200m and 
the mass enclosed in it as M2oom- We consider all the halos hav- 
ing A/2oom > 10 14 Mq h~ x . In addition, we selected subsam- 
ples of the 200 objects with Af 2 oom closest to 7 x 10 13 M hT 1 , 
5 x 10 13 M Q h~\ 3 x 10 13 M© ft -1 , and 10 13 Mq h~\ Starting 
from the centres of the halos, we construct radial profiles by bin- 
ning the particles in radial bins. We select and study objects at three 
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Figure 2. The values of x a f f in units of i?200m for the objects in ACDM 
at 2 = 0. The vertical dashed blue line corresponds to the value defining 
relaxed objects, x tf = 0.07-f?200m- The vertical pink line marks the 



median value, x 



off = 0.046i?200r, 



Figure 3. The values of oy ma for the objects in ACDM at z 
vertical pink line marks the median value, a rms = 0.0358. 
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different redshifts: z — 0, z = 0.5, and 2 = 1. For the following 
analysis, we also calculate for each cluster selected in this way the 
radius at which the overdensity drops below 200 (500) times the 
critical background density and denote it as ifeoo (Rsoo)- The cor- 
responding mass is indicated as A/200 (-Msoo). It is useful to define 
a quantitative criterion to decide whether a cluster can be consid- 
ered relaxed or not because, in general, relaxed clusters have more 
spherical shapes, better defined centres and thus are more repre- 
sentative of the self-similar behaviour of the dark matte r halos. We 
use a simple criterion similar to the one introduced in iNeto et al.l 
(2007): first we define x f / as the distance between the centre of 
the halo (given by the most bound particle) and the barycentre of 
the region included in i?200m ; then we define as relaxed the halos 
for which x // < 0.07i?200m. We plot the distribution of x ff for 
the objects in ACDM at 2 = in Fig.[2] 



3 FITTING c M RELATION 

In this paper, we analyse both the dark matter only and the hy- 
drodynamical runs of the simulation set introduced in Paper I. We 
calibrate the fitting procedure on ACDM because we use it as a 
reference model, and most references in the literature are based on 
this model. 

For the concordance ACDM model, for each cluster at 2 — 
in the dark matter only run, we perform a logarithmic fit, using 
Poissonian errors (In 10 x y/ndm) -1 (where ridm is the number 
of dark matter particles in each radial bin, of the order of 10 — 10 3 
depending on the mass of the object), of the three-dimensional dark 
matter profile pdm(r) in the region [0.1 — l]i?200 (where the value 
of R20 is taken directly fr om the true mass profile) with a NFW 
profile jNavarro et aTll 19961) 



S = 



200 



c 200 



pdm (r) 



(r/r s )(l + r/r s ) 2 
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where p c is the critical density, r s is the scale radius and 5 is a 
characteristic density contrast. Then, instead of defining C200 = 
R20o/r s , we directly find the concentration parameter C200 from 
the normalization of the NFW profile 



[ln(l + C200) 



1 + C200 
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We require the central density parameter 8 to be greater than 100 
and the scale radius r s to be within [0.1 — ljifeoo- We exclude the 
inner regions from the fit because we are limited in resolution inside 
a given radius. We indicate the dark matter concentration found in 
this way as C2oodm- We define the rms deviation oy ms as 



N bins 

■ ^2 [lOgioPi - logiopjVFw]" 
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where Nbins is the number of radial bins over which the fit is per- 
formed. We plot the distribution of a rms for the objects in ACDM 
at 2 = in Fig. [3] 

In addition to the complete and relaxed samples discussed in 
Sect. [2] we create a sample of "super-relaxed" objects defined as 
the ones having both x ff and ay ms lower than the median value 
extracted from the complete sample, namely 0.046 and 0.0358, re- 
spectively. We do this to check if strict restrictions on both the dy- 
namical state and the shape of the profile of the objects can reduce 
the intrinsic scatter in the values of concentration. We stress that the 
definition of super-relaxed objects is dependent on the way we fit 
the profile. In the end, the complete sample, the relaxed sample and 
the super-relaxed sample at 2 = are constituted by 1357, 923, 
and 41 1 objects, respectively. 

We bin the objects in the complete sample in groups of 200, so 
that we have bins around 10 13 M ft -1 , 3 x 10 13 M Q h' 1 , 5 x 

10 13 Mq /i-\and7x 10 13 M Q ft" 1 . For halos more massive than 

10 14 Mq hT 1 , we bin the objects starting from the low-mass ones, 
so that the most massive bin can contain less than 200 objects. The 
analysis for the relaxed and super-relaxed samples is done selecting 
the relaxed and super-relaxed objects inside each bin. Once we have 
C200dm for each object in each mass bin, since the distribution of 
C200dm is log-normal inside each bin, we evaluate the mean M200 
and the mean and rms deviation of logioC2oodm in each bin, for all 
the three samples. In the following of the paper, when we indicate 
the value of C2oodm in a mass bin, we refer to l0^ logloC200dm \ 

In Fig. [4] we plot C2oodm for each object in the complete sample. 
We see that there is a large intrinsic dispersion in the values of the 
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Figure 4. The values of c 2 oodm f° r me complete sample of the ACDM 
model at z = 0. For each object, we plot M200 an d c 200dm- Differ- 
ent colours indicate different M200m ranges. The black squares indicate 
c 200dm m eac h mass bin (see text for details). For each mass bin, we plot 
the mean M200 and C2oudm with rm s deviation. 



Table 2. Best-fit parameters, standard errors and reduced chi-squared \ 2 
of the c — M relation equation 4 1 01 for dark matter only density profile fit 
in the region [0.1 — l]i?200 f° r the complete, relaxed and super-relaxed 
samples of the ACDM model at z = 0. 



Model 


C8 


A 




B 




x 2 


ACDM 


0.776 






dm 






all 

relaxed 
super-relaxed 




3.59 
4.09 
4.52 


0.05 
0.05 
0.06 


-0.099 
-0.092 
-0.091 


0.011 
0.011 
0.013 


0.48 
0.66 
0.76 



mass. Since we want to compare the dark matter only runs with the 
hydrodynamical ones, and in the hydrodynamical runs we do not 
completely resolve the baryonic physics on very small scales, we 
will take a conservative approach and fit in the range [0.1 — l]i?200 • 
Still, it is important to know what happens if we consider also the 
inner regions, in particular if we want to compare our results with 
the ones in literature. 

With the mean and rms deviation of logioCsoodm in each bin at 
hand, we fit, for the complete, relaxed and super-relaxed samples, 
the binned c — M relation using 



concentration inside each group of objects, which is marked with a 
different colour. For objects more massive than 10 14 M h' 1 , the 
maximum value of c 2 oo<im clearly decreases with M200 • Only when 
we bin the objects and plot the mean value of C2oodm in each mass 
bin we can see a general trend with the concentration decreasing 
with increasing mass, even if the rms deviations are quite large. 
In Fig.|5]we plot the values ofc2oodm in each mass bin for the com- 
plete, relaxed and super-relaxed samples, along with the number of 
objects in each bin. We note that, in all bins, more than 50% of the 
halos are relaxed and around 30% are super-relaxed. In each bin, 
the value of C2oodm for the relaxed sample is higher than the one 
for the complete sample, and the value for the super-relaxed sample 
is even higher. In all three samples, C2oodm is decreasing with in- 
creasing mass. In Fig. [5] we also show the scatter in logioC2oodm in 
each bin for the complete, relaxed and super-relaxed samples. For 
the complete sample, the scatter ranges from 30% to 40%, with a 
positive trend with mass. It means that inside each mass bin there 
are objects with quite different concentrations, in particular in the 
high-mass tail. For the relaxed sample, the scatter reduces to 20% 
up to 30%, meaning that part of the scatter in the complete sam- 
ple is due to objects in a particular dynamical state. If we move to 
the super-relaxed sample, we notice that the scatter stays between 
15% and 20%. On the one hand, this means that putting strong con- 
straints on both the dynamical state and the shape of the profile of 
the halos allows us to halve the intrinsic scatter; on the other hand, 
even when considering the most relaxed and smooth objects in our 
sample we cannot reduce the intrinsic scatter below 15%. 
As a check, for the complete sample of the ACDM model at 2 = 
we also evaluate c 2 oodm by fitting the NFW profile equation 10 
in the range [0.01 — ljifeoo- For low mass objects in our sample, 
the value of O.OI-R200 is close to the limit of the force resolution 
of the simulation, which is 7.5 kpc h^ 1 at z = 0. In Fig [6] we 
show the results for the values of C2oodm and the relative error in 
logioC2oodm in each mass bin. We see that, by fitting in the range 
[0.01 — 1]7?200, we obtain concentrations up to 10% higher than 
by fitting in the [0.1 — l]i?200 range. Moreover, by fitting includ- 
ing the inner regions in the fit, the intrinsic scatter in concentra- 
tion is lower by about 5%. These trends are almost independent of 



logioc 2 oo = logioA + B logio 



10 14 M 0j 



(10) 



where logioC2oo and M200 are the mean values in each bin. For the 
error on the mean of logioC2oodm in each bin, o"c, we use the rms 
deviation of logioC2oodm divided by the square root of the number 
of objects in the bin. For each fit we also define 
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logl C200i — logl0C 2 00 /it 
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and evaluate the reduced chi-squared x 2 , i.e. x 2 divided by the 
number of degrees of freedom. We list the best fit values A and 
B for each sample, along with the corresponding standard errors 
and \ 2 , in Tablef2] We see that, compared to the complete sample, 
the normalization A increases by about 15% for the relaxed sample 
and by about 25% for the super-relaxed sample, while the slope B 
does not change significantly, even if excluding unrelaxed objects 
results in a shallower slope. The values of \ 2 indicate that equation 
d 1 0b is a good parametrization of the c — M relation in logarithmic 
scale. 

In order to understand the impact of low-mass object on the c — M 
relation, we check how the best-fit values of the c — M relation 
change if we do not include the less massive objects. We report 
the results we obtain by considering only objects with M200W1 > 
10 14 Mq ft -1 in Table |3] For all three samples, we find a flatter 
relation than when including also low-mass objects, with larger er- 
rors on the slope. The normalizations are lower of few percentage 
points, while the relative errors are a factor of two higher compared 
to the case where low-mass objects are also considered. Moreover, 
we find that in this case the slope is very sensible to the dynamical 
state of the objects included in the sample. Thus we can conclude 
that the inclusion of low-mass objects is necessary to find a sig- 
nificant correlation between the concentration and the mass of the 
halos. We do not quote the reduced chi-squared in this case because 
just 3 mass bins are considered with 2 parameters to be fitted. 

For the dark matter only profiles of the concordance ACDM 
model we also perform a logarithmic fit of equation (|7) without 



6 C. De Boni et al. 



6.0 
5.5 
5.0 
E 4.5 

■o 

O 

o 
fN 

° 4.0 

3.5 

3.0 
2.5 



□ 

40 

□ 



200 
□ 



10' 



ACDM z=0 



■-■i 
□ 



1 37 
□ 

200 
□ 



oil 

relaxed 
super-relaxed 



□ g 

138 129 

□ □ 



□ 



□ 



200 
□ 



000 

□ 



133 13S 
□ □ 



200 
□ 



ooo 
□ 



M 200 ( M o) 



10" 



□ 



I 00 

□ 



□ 



ACDM z=0 



0.4 



o 0.3 



0.2 



0.1 



• oil 

• relaxed 

• super-relaxed 



137 138 129 



10' 3 10"* 

M 200 ( M o) 
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Table 3. Best-fit parameters and standard errors of the c — M relation equa- 
tion (10) for dark matter only density profile fit in the region [0.1 — 11-R.200, 
considering only objects with M^oOm > 10 14 Mq h~ 1 , for the complete, 
relaxed and super-relaxed samples of the ACDM model at z = 0. 



Model 




A 




B 




ACDM 


0.776 






dm 




all 

relaxed 
super-relaxed 




3.55 
3.99 
4.35 


0.09 
0.10 
0.13 


-0.087 
-0.055 
-0.010 


0.038 
0.042 
0.049 



using Poissonian errors, as usually found in the literature. In this 
case we evaluate C200 both from equation <[8l» and by directly defin- 
ing C200 = R20o/r s (using -R200 from the true mass profile), and 
indicate the two values as C2oodm,fit and C2oodm,rec, respectively. 



Moreover, in order to check the robustness of our fit, we perform 
a logarithmic fit, using Poissonian errors, of the dark matter profile 
times r 2 , and a logarithmic fit, using Poissonian errors, of the dark 
matter profile times r 3 . 

Finally, we perform a logarithmic one-parameter fit, using Poisso- 
nian errors, of equation (|7} re-expressed as 

Pc 

= 200 cjoo 1 

3 [Ml + «*»)-!&] fe][l+fe)] 2 ' 

where the only free parameter is C200, since R200 is directly taken 
from the true mass profile. 

Only f or comparison, we also evaluate the halos concentration fol- 
lowing IPradaetaljHoT 1), i.e. by solving 



c — M relation for galaxy clusters 



ACDM all z = 



where 



V200 



Vmax — max 



0.216 c 

~1W 



1/2 



GM{< r) 



1 1/2 



t r ( GM200 \ 

1/200 = 

and, for a NFW profile, /(c) is given by 



/(c) = In(l + c) 



1/2 



1 + c 



(13) 



(14) 



(15) 



(16) 



We find that the best values [8, r s ], and therefore C200 from equa- 
tion (8), found by fitting the dark matter profile times r 2 and the 
dark matter profile times r 3 are exactly the ones found from equa- 
tion (|7) with Poissonian errors, while the ones from equation {7J 
without Poissonian errors, from equation ( 1121 . and from equation 
d 1 3 fc are somehow different. We show these differences for the com- 
plete, relaxed, and super-relaxed samples in Fig. [7j For the com- 
plete sample, we notice from Fig. [7] that the unweighted fit tends 
to give values of C200dm slightly higher compared to the Poisson 



weighted fits, with C2oodn 



always higher than c 2 ood 



m,fit 



The 



values of 0200dm found from the one-parameter fit of equation (112) 
are consistent with the ones found from equation (|7} with Poisso- 
nian errors. All these models are compatible within f ew percent- 
age po ints, while the method using Vmax discussed in lPrada et al.l 
J201 II) gives systematically higher values, with a positive trend 
with mass. In the most massive bin, the difference between the 
concentrations is more than 20%. For the relaxed sample, the sit- 
uation is similar, but the differences among the different fits are 
somewhat smaller. Even for the method based on equation dl3h 
which is the one that gives very diverse results, the difference is 
at most around 10%. For the super-relaxed sample, the differences 
among the various fits almost disappear, with the exception of the 
method using Vmax, which however now shows differences below 
5% and, in the lowest mass bin, recovers values of the concentra- 
tion even lower than the other methods. The method based on Vmax 
recovers the same values of the concentration of the usual NFW fit 
for super-relaxed objects, as expected because it is based on the 
NFW profile. But for less relaxed objects it tends to maximize the 
value of Vmax /V200 by taking into account features due to accre- 
tion/merging, thus artificially enhancing the concentration of the 
halos. 

In Table|4]we show the best-fit parameters, the standard errors 
and the reduced chi-squared of the c — M relation for the different 
fits discussed above. As expected, for all the three samples, the dif- 
ferences in the normalization and slope are limi ted to a few precen t, 
with the exception of the method discussed in lPrada et al.l (201 1). 
For this method we find a flatter relation, due to the fact that it 
predicts higher concentrations for high-mass objects compared to 
the usual fit of the NFW profile equation (0. Moreover, also the 
normalization is higher (more than 15% for the complete sample, 
more than 5% for the relaxed sample, and less than 5% for the 
super-relaxed one), given that the concentration found from equa- 
tion l |13t is in general higher than the one found in the usual way. As 
already shown, the effect is stronger for the complete sample than 
for the relaxed and super-relaxed ones. The values of the reduced 
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Figure 7. (Top panel) Ratio between c 2 oodm (black squares), c 2 oodm,fit 
(blue triangles), C2oodm,rec (blue diamonds), C2oodm found by fitting 
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Table 4. Best-fit parameters, standard errors and reduced chi-squared \ 2 of 
the c — M relation equation HOI for the different fits of the dark matter only 
density profile for the complete, relaxed and super-relaxed samples of the 
ACDM model at z = 0. 



ACDM z=0 



Model 


0"8 


A 




B 


° B 


A 


ACDM 


0.776 






dm (all) 






Poisson 




3.59 


0.05 


-0.099 


0.011 


0.48 


c 200dm,fit 




3.69 


0.04 


-0.102 


0.010 


0.49 


^200dm,rec 




3.71 


0.04 


-0.101 


0.010 


0.57 


one-parameter 




3.61 


0.05 


-0.096 


0.011 


0.50 






4.18 


0.04 


-0.063 


0.008 


0.67 


ACDM 


0.776 






dm (relaxed) 






Poisson 




4.09 


0.05 


-0.092 


0.011 


0.66 


C200dm,fit 




4.14 


0.05 


-0.096 


0.010 


0.53 


^200dm,rec 




4.17 


0.05 


-0.095 


0.010 


0.55 


one-parameter 




4.12 


0.05 


-0.088 


0.011 


0.74 






4.40 


0.04 


-0.065 


0.008 


0.90 


ACDM 


0.776 




dm (super-relaxed) 




Poisson 




4.52 


0.06 


-0.091 


0.013 


0.76 


C20Qdm,fit 




4.55 


0.06 


-0.092 


0.013 


0.67 


^200dm,rec 




4.56 


0.06 


-0.091 


0.013 


0.66 


one-parameter 




4.56 


0.06 


-0.086 


0.012 


0.63 


Vmax 




4.69 


0.06 


-0.060 


0.010 


0.45 



chi-squared indicate that equation dlOt is a good parametrization of 
the c — M relation also for these different definitions of the con- 
centration. 

Before moving to the hydrodynamical simulation, we sum up here 
our findings for the dark matter only run of the reference ACDM 
model. The intrinsic dispersion in the logarithmic values of C200dm 
for objects of similar mass is reduced by a factor of two if we limit 
our analysis to objects that are both relaxed and with a clear NFW- 
like shape of the dark matter profile. Still, at best, the intrinsic scat- 
ter is of the order of 15%. The more relaxed the objects in the sam- 
ple, the higher the normalization A of the c — M relation, while 
the slope B is almost independent of the dynamical state of the ha- 
los. This is no longer true if we focus our analysis to objects with 
M2oom > 10 14 Mq hT 1 , for which the slope is shallower for the 
relaxed sample and almost flat for the super-relaxed one. If we limit 
ourselves in fitting the dark matter profile with a NFW profile, we 
almost recover the same values of C2oodm , A, and B, independently 
of the way we treat the errors on the fit and on the number of free 
parameter s that we fit . Thing s do change if we use the method dis- 
cussed in lPrada et al.l d201ll) . which systematically overestimates 
the concentration compared to the others, in particular for high- 
mass objects. This results in both higher values of A and B for this 
method. With this in mind, in the following of this paper, in partic- 
ular when we study the effect of dark energy on the c — M relation, 
we will always consider Cwodm recovered from equation l(8) from 
the fit of equation Q with Poissonian errors, and distinguish only 
between the complete and the relaxed sample, whose definition is 
independent of the way we find the concentration. 



3.1 Hydrodynamical runs 

For the hydrodynamical run, we fit the total three-dimensional den- 
sity profile ptot = Pdm + Pgas + p star in the range [0.1 - lJ-Raoo 
with equation {7}, using Poissonian errors (In 10 x \frii) ~ 1 for each 
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Figure 8. Ratio between c 2 oodm (squares) and C2ootot (stars) for the com- 
plete (black), relaxed (blue) and super-relaxed (pink) samples and C2oodm 
for the complete sample of the ACDM model at z = 0. 



component (where m is the number of particles of the i-th species) 
and summing them in quadrature. We indicate the total matter con- 
centration found from equation ([8} with C2ootot- We show the dif- 
ferences between C2oodm and C2ootot for the complete, relaxed and 
super-relaxed samples in Fig. [8] We clearly see that, in all three 
samples, starting from objects with M^oo-m > 10 14 Mq the 
concentration in the hydrodynamical run is higher than in the dark 
matter only case, and the effect becomes more relevant at higher 
masses. This effect is less pronounced in objects with M^oom < 
10 14 Mq ft -1 . Thus, the inclusion of baryons appears to affect 
more massive galaxy clusters than small groups. In order to check 
this fact, we evaluate the relative distribution of baryons and stars 
inside -R200 and in the range [0.1 — 0.3]-R2no, i-e. the innermost part 
of the range in which we fit the NFW profile. We show the results 
in Fig. [9] We see that, while the total baryon fraction at -R200 is 
almost constant with mass, if we limit to the range [0.1 — 0.3]i?2oo 
the baryon contribution to the total mass becomes more important 
for more massive objects. The same happens for the stars, which 
are known to concentrate in the internal regions of halos. Thus, in 
this simulation, the relative contribution of baryons and stars in the 
inner regions is more relevant in massive galaxy clusters than in 
small groups, and this fact leads to an increase of the concentration 
in massive objects when including baryonic physics in the simula- 
tions. 

Then, we fit the c — M relation equation dlO) as in the dark matter 
only case. We list the best-fit parameters, standard errors and re- 
duced chi-squared of the c — M relation for the complete, relaxed, 
and super-relaxed samples in Table [5] For all samples, the normal- 
ization A is from 5% to 10% higher compared to the dark matter 
only case, while the slope is shallower by about 30%. This is ex- 
pected because, as we have already seen, in the hydrodynamical 
runs high-mass objects are more concentrated than objects of the 
same mass in the dark matter only simulations. For the hydrody- 
namical runs the values of the reduced chi-squared are quite high. 
This fact seems to indicate that baryons introduce some dependence 
and equation <10K even though it remains a good parametrization, 
is no more able to completely characterize. 
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Figure 9. (Top panel) Star fraction fstar (red) and total baryon fraction 
fbar (black) at -R200 (points) and in the [0.1 — 0.3]i?,2OO range (triangles) 
for the complete sample of the ACDM model at z = 0. (Bottom panel) 
Ratio between f s t ar in the range [0.1 — 0.3]iJ200 and f s tar at -R200 (red 
triangles) and between fi, ar in the range [0.1 — 0.3]-f?200 and/f, ar at_R200 
(black triangles) for the complete sample of the ACDM model at z = 0. 



Table 5. Best-fit parameters, standard errors and reduced chi-squared x 2 
of the c — M relation equation j 101 for dark matter only and total density 
profile fit in the region [0.1 — l]-f?200 f° r tne complete, relaxed and super- 
relaxed samples of the ACDM model at 2 = 0. 



Model 


0"8 


A 




B 




x 2 


ACDM 


0.776 






all 






dm 
total 




3.59 
3.81 


0.05 
0.05 


-0.099 
-0.061 


0.011 
0.011 


0.48 
1.69 


ACDM 


0.776 






relaxed 






dm 
total 




4.09 
4.29 


0.05 
0.05 


-0.092 
-0.064 


0.011 
0.011 


0.66 
1.54 


ACDM 


0.776 






super-relaxed 






dm 
total 




4.52 
4.89 


0.06 
0.06 


-0.091 
-0.062 


0.013 
0.011 


0.76 
1.34 
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Figure 10. Ratio between c 2 oodm (squares) and C200tot (stars) for the com- 
plete (black) and relaxed (blue) samples and C2oodm f° r the complete sam- 
ple of the ACDM model at z = 1. 



3.2 The dependence upon redshift 

Finally, in order to study the evolution with redshift of the c — M 
relation, which is of fundamental importance if we want to distin- 
guish among different cosmological models, we also consider ob- 
jects at z — 0.5 and z — 1. We show the differences between 
C200dm and C2ootot for the complete and relaxed sample of the 
ACDM model at 2 = 1 in Fig. \l0\ Even if the trend with mass 
is less clear than at 2 = 0, still we can see that, already at z = 1, 
C200tot is greater than c 2 oodm in both samples. 
When we consider also objects at 2 = 0.5 and z = 1, we fit a gen- 
eralized form of the c — M relation equation dlOt with an explicit 
redshift dependence, namely 



logl0C200 = 

logicM + B logio 



(17) 



M 2 . 



10 14 M 



Clogi (l + 2) 



We can perform this fit in two ways. Either we keep the best-fit 
values Aq and Bo found at 2 = fixed and fit only C in equation 
dl7K or we perform a three-parameter fit by keeping A, B and C 
free. We report the results of the different fits both for the complete 
and relaxed samples of the dark matter only and hydrodynamical 
simulations in Table [6] We note that, in all cases, the redshift de- 
pendence is negative, meaning that objects of a given mass have 
lower concentration at higher redshift. This is expected, because 

— 1/3 

c ~ p c and the critical density drops with decreasing redshift. 
If we leave the normalization and the slope free, we see that A 
changes at most of few percentage points, while B can vary sig- 
nificantly. The redshift dependence seems to be insensitive to both 
the way in which the fit is performed and the dynamical state of 
the halos, while it is different for dark matter only and total con- 
centration, being steeper in the former case and shallower in the 
latter. In general, the reduced chi-squared of the fit of equation dl7| ) 
is rather high, in particular for the hydrodynamical runs. This can 
be an indication that the redshift dependence we are considering is 
somehow too simple to fully reproduce the redshift evolution of the 
c — M relation, and that the presence of baryons makes this evo- 
lution more complex. In the following sections, we will constrain 
the normalization and slope at 2 = 0, and then study the redshift 
evolution keeping A and B fixed. 
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Table 6. Best-fit parameters, standard errors and reduced chi-squared \ 2 of the c — M relation equation Q 71 for dark matter only and total density profile fit 
in the region [0.1 — 1J-R200 f° r the complete and relaxed samples of the ACDM model at z = 0, z = 0.5 and z = 1 both fixing A and B at the best-fit 
values at z = and keeping all the parameters free. 



Model 


CT8 


A 




B u B 


C 




x 2 


ACDM 


0.776 






dm (all) 








fixed 
free 




3.59 
3.63 


0.04 


-0.099 — 
-0.077 0.008 


-0.33 
-0.32 


0.02 
0.03 


1.68 
1.41 


ACDM 


0.776 






dm (relaxed) 








fixed 
free 




4.09 
4.13 


0.05 


-0.092 — 
-0.080 0.008 


-0.31 
-0.31 


0.02 
0.03 


1.22 
1.24 


ACDM 


0.776 






total (all) 








fixed 
free 




3.81 
3.81 


0.04 


-0.061 — 
-0.046 0.007 


-0.26 
-0.25 


0.02 
0.03 


2.40 
2.40 


ACDM 


0.776 






total (relaxed) 








fixed 
free 




4.29 
4.27 


0.05 


-0.064 — 
-0.071 0.007 


-0.26 
-0.26 


0.02 
0.03 


1.97 
2.15 



4 MODEL-INDEPENDENT CONCENTRATION 

All the fits we used in the previous sections are based on some 
assumptions that we made on the shape of the density profile. In 
principle, it should be useful to have a model-independent proxy 
for the concentration found from equation (|8}. Given the usual def- 
inition of the concentration, a natural quantity could be the ratio 
between two physical ra dii, characterized b y a given overdensity. 
As already pointed out in lDuffv et al.|j201ol) . a good choice should 
be i?2500 / ^?500 since these radii are commonly used in X-ray ob- 
servations. For each object, we evaluate the ratio -R2500 /Rmo tak- 
ing -R2500 and Rmo directly from the true mass profile, without 
any assumption on the density profile. We bin the objects in mass 
as in the previous sections. The results for the complete, relaxed 
and super-relaxed samples, are shown in Fig. Qj] We note is that, 
for a given sample, the variation of C200 is more than 30% while 
the variation of -R2500 /R500 is less than 10%. This means that the 
ratio between -R2500 and -R500 is less dependent on the mass of 
the halo when compared to C2oo- We also notice that objects with 
higher concentration, i.e. with higher R20o/ r a , also have higher 
i?250o/-Rsoo. This result confirms the ones in iDuffv et alj (2008), 
where a weak positive trend of -R500 /-R2500 with M500 was found. 
Finally, there is a clear dichotomy between the dark matter only 
and the total profile fits. At a given concentration, when including 
baryons, the ratio -R2500AR500 is more than 5% higher with com- 
pared to the dark matter only case. Moreover, in the dark matter 
only case, the relation between R25001 'Rsoo and C200 closely fol- 
lows the one predicted by a NFW profile, indicating that the halos 
are well described by this profile. The different behaviour of the 
halos extracted from the hydrodynamical run suggests that, when 
we consider the total (dark matter plus baryons) density profile, the 
NFW profile is no longer a good approximation to the real profile. 



5 COMPARISON WITH OTHER WORKS 

In this section, we compare with previous works our results on the 
dark matter concentration obtained by fitting equation in the 
range [0.1 — ljifeoo- We use objects at z = 0, z = 0.5 and 
z — 1 from the dark matter only simulations of the concordance 
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Figure 11. Comparison between C2oodm (squares) and C2ootot (stars) and 
^50o/^?2500 f° r the complete (black), relaxed (blue) and super-relaxed 
(pink) samples of the ACDM model at z = 0. The black dashed-line repre- 
sents the prediction for a NFW profile. 



ACDM model normalized with WMAP3 data. We recall that, by 
fitting the c — M relation equation d!7t . we find that more mas- 
sive objects are less concentrated than less massive ones and that 
objects at high redshift are less concentrated than objects at z = 0. 
Moreover, relaxed objects are more concentrated compared to the 
complete sample. All these findings qualitatively confirm what can 
be found in literature. For the sake of comparison with literature, 
some numbers are summarized in Table [7] where we have adjusted 
the quoted normalization A to our pivot mass 10 14 Mq using the 
correspondin g quoted slope B. 

Mac cio et alj(2008) make a comparison of the c — M relation for 
all and relaxed objects in different ACDM cosmologies, namely 
the ones obtained using the parameters coming from WMAP1, 
WMAP3 and WMAP5. They fit the dark matter profile with a 
NFW profile. They span the mass range 10 10 M h' 1 < M < 
10 Mq ft -1 and bin the objects in mass bins of 0.4 dex width. 
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Table 7. Best-fit parameters of the c — M relation, comparison with other works. See text for details. 



Reference 


h 




C8 




dm (all) 






dm (relaxed) 










A 


B 


C 


A 


B 


C 


This work (fixed A and B) 


0.704 


0.268 


0.776 


3.59 


-0.099 


-0.33 


4.09 


-0.092 


-0.31 


This work (free A and B) 


0.704 


0.268 


0.776 


3.63 


-0.077 


-0.32 


4.13 


-0.080 


-0.31 


Maccio et al. (2008) 


0.71 


0.268 


0.90 


4.55 


-0.119 




5.31 


-0.104 






0.73 


0.238 


0.75 


3.60 


-0.088 




4.13 


-0.083 






0.72 


0.258 


0.796 


3.84 


-0.110 




4.46 


-0.098 




Neto et al. (2007) 


0.73 


0.25 


0.9 


4.85 


-0.11 




5.45 


-0.10 




Duffv et al. (2008) 


0.742 


0.258 


0.796 


4.06 


-0.097 




4.81 


-0.092 




0.742 


0.258 


0.796 


4.23 


-0.084 


-0.47 


4.85 


-0.091 


-0.44 


Pradaetal. (2011) 


0.70 


0.27 


0.82 


5.31 


-0.074 




5.55 


-0.08 





WMAP3 and WMAP5 data suggest a lower matter density Q m and 
a lower power-spectrum normalization erg than WMAP1 ones. This 
means that halos of a given mass form later, and thus should be less 
concentrated. And indeed this is what they find, as it can be seen 
from Table [7] The higher <rg and fi m , the higher the normalization 
and the steeper the slope of the c—M relation. With respect to their 
results for the WMAP3 cosmology, we find the same values of the 
normalization both for the complete and relaxed samples, with a 
slightly steeper slope. We show this excellent agreement in Fig.l 121 

In a more recent paper, iMunoz-Cuartas etaL I fcOllh fit the c- M 
relation using a formula similar to equation ( 1 101 . but incorporating 
the redshift dependence by letting A and B to be functions of red- 
shift themselves. For the c v j r — M v i r relation, a t z = they are able 
to reproduce the results of lMaccio et al.l d2008l) with differences of 
the order of few percentage points. 

iNeto et al.1 ( [200% make a compari son, using halos extra cted from 
the Millennium Simulation (MS) ISpringel et al. 2005), between 
the complete sample and relaxed objects only. They fit the dark 
matter profile with a NFW profile in the [0.05 — l]R V ir range, with- 
out Poisson weighting. They span the mass range 10 12 Mq h^ 1 < 
M < 10 15 Mq h^ 1 and bin the objects in mass. They find that re- 
laxed objects have more than 10% higher normalization compared 
to the whole sample, while the slope is 10% shallower. We find 
the same trend, but there is a discrepancy in the absolute numbers, 
both in the complete and relaxed samples. Indeed, they find nor- 
malizations that are more than 30% higher than ours, which can 
be explained with the higher <rg (but also with the different radial 
range over which they fit, see Fig.|6]in Sect. [3}, and a 10% steeper 
slopes. 

iDuffv et all J2008h make a comparison with lNeto et all d2007l) us- 
ing the parameters suggested by WMAP5 (CMB only) data. They 
fit the dark matter profile with both a NFW and an Einasto pro- 
file in the [0.05 — l]R V ir range. We discuss the results for the 
NFW profile case. They span the mass range 10 11 Mq ft -1 < 
M < 10 15 Mq h -1 a nd bin the objects in mass. In comparison 
with INeto et alj (2007), they find lower normalizations by about 
15% due to a lower value of erg. They also fit the c — M relation 
taking into account the redshift evolution from z = to z — 2, 
both for the complete sample and the relaxed objects only. They 
find a stronger dependence on redshift than what we find, both for 
the whole sample and relaxed systems. A possible explanation is 
that we fit in the range between z — and 2 = 1, while they 
reach z = 2. If the redshift evolution is not constant with red- 
shift, but it is weaker at low redshift, this could be a possib le ex- 
planation for the different slope we find. lDuffvetai] (2008) also 
find that the concentration of their halos is lower than the one in- 



ferred from X-ray observations and ascribed that fact to the effect 
o f baryon physics th at was missing in their simulations. However, 
in lDuffvetal.l ( f2010l) it is shown that even including baryon physics 
in the simulations they cannot reproduce both observed concentra- 
tions and stellar fraction in galaxy groups and clusters. In general, 
at z = 0, when including metals and AGN feedback, they find 
lower concentrations compared to the dark matter only case, and 
the effect is more relevant in low-mass objects. When including 
only primordial cooling, they find higher concentrations compared 
to the dark matter only case, qualitatively in agreement with what 

we find in F ig. [8] 

[ Prada et all d201 lh . using the Bolshoi simulation ( Kly pin et al.l 
1201 lh . evaluate the concentration of the halos using equation J13l >. 
binning the objects in V ma x and fitting a c — <j(M, z) relation, 
where <j(M, z) is the linear rms fluctuation of density field on the 
scale M. They find that the concentration c(<j) has a nearly univer- 
sal U-shaped profile, with some small dependence on redshift and 
cosmology. They also provide a fit of equation dlOt for all their ha- 
los at z = 0. They find a higher normalization and a shallower slope 
in comparison with other works, as we do when we use equation 
dl3K instead of equation {8}, to evaluate the concentration. For re- 
laxed halos, selected by V ma x , they find a 5% higher normalization, 
as we do. Moreover, they find that the differences in concentrations 
for relaxed halos selected by V ma x are higher for high-mass objects 
than for low-mass ones. 

We conclude this section noting that, when the values of the 
cosmological parameters are similar, our findings about the c — M 
relation in the reference ACDM model are in good agreement with 
what already found in literature. So we can safely rely on our 
ACDM model as a reference, when comparing the impact of dif- 
ferent dark energy models on the c — M relation. 



6 DARK ENERGY MODELS: RESULTS ON THE DARK 
MATTER PROFILES 

From now on, we start to compare the c — M relation in the ACDM 
cosmology with the ones derived for the other cosmological models 
under investigation. The c—M relation for galaxy clusters extracted 
from dark matter only simulations of different dark energy models , 
including RP and SUGRA, has been studied in lDolag et ail b004l) . 
They fit a formula similar to equation ( 117) , but keeping C fixed to 
— 1. They find that, when the same <rg is used for all the models, 
the c — M relations for dark energy cosmologies have higher nor- 
malizations compared to ACDM, depending on the ratio between 
the growth factors through 
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ACDM all z = ACDM relaxed z = 
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Figure 12. (Left panel) The values of c 2 oo£j m for the complete sample of the ACDM cosmology at z = 0. The black dotted-dashed line is our best-fit of 
c — M relation equation (To) and the vertical black bar is the error on the normalization as listed in Table|7] The pink solid line is the best-fit relation of Maccio 
et al. (2008) for the WMAP3 cosmology. (Right panel) The same as in the left panel, but for the relaxed sample. 
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where the co llapse red s hifts z co u are evaluated following the pre- 
scriptions of Ek e et al.l J200lh . When as values are normalized to 
CMB data, as we do in this work, dark energy cosmologies have 
lower normalizations compared to ACDM. We find that, in order 
to recover the values of the normalization they quote in this case, 
equation j I Sl i should be multiplied by the ratio between the values 
of as, i.e. cts. de /cs.acdm- This f act goes in the same direction as 
what found in lMaccio et al.l ( 120081) , where models with higher as 
also have a higher normalization of the c — Al relation. 

We recall that to obtain the concentration we fit equation 
in the range [0.1 — l]i?200 using Poissonian errors and adopting the 
best-fit parameters to obtain C200 from equation l[8}. Then we bin 
the objects in groups of 200 halos, starting from less massive ones, 
and also define a relaxed sample by taking the relaxed objects inside 
each bin. Finally, we fit the binned c — M relation with equation 
j 101 >. We begin the comparison in this section with the dark matter 
only runs at z = 0, while in the following section we will study the 
hydrodynamical runs, also at higher redshifts. 

In Table [8] we summarize the best-fit parameters, the standard 
errors and the reduced chi-squared of the c — M relation equa- 
tion J 10b for the five cosmological models here considered, both 
for the complete and relaxed samples. For the complete sample, the 
differences in the normalization A between ACDM and the other 
cosmological models are less than 10%, with EQn being the only 
model having a higher normalization. The slope B is within 5% 
of the ACDM value for all the models with the exception of EQn 
again, which shows a 30% flatter slope. For the slope the differ- 
ences among the models, excluding EQn, are smaller than the stan- 
dard errors, while for the normalization these differences are sig- 
nificant. If we limit ourselves to the best-fit values, given that the 
slopes are almost identical and that all the cosmological parameters 
except (78 are fixed, we expect that the normalizations should fol- 
low the values of as, i.e . the higher as the higher the normalization 
(see iMaccio et al]|2008l) . an d D+, i.e. the high er D+ at z co u the 
higher the normalization (see D olag et al.l2004l) . The quantity con- 
trolling the normalization is thus expected to be asD+(z co u). We 



plot the ratio between the value of asD+ for a given dark energy 
model and the one for ACDM as a function of redshift in Fig. [T3l 
Independently of the precise definition of z co u, the cosmological 
model with the highest value of this quantity is ACDM, followed 
by RP, EQp, EQn, and SUGRA. We do expect the normalizations 
of the c — M relation to follow the same order, with ACDM hav- 
ing the highest and SUGRA the lowest. Instead we see that, on the 
one hand, EQp which has the third highest asD+ has the lowest 
normalization while, on the other hand, EQn which has the sec- 
ond lowest (78-D+ has the highest normalization. The relative order 
of asD+ and A i s preserved for ACDM, RP and SUGRA, as in 
iDolag et al.l j2004) . 

In a recent paper, Bhattacha rva et all d201 lh . using Af-body numer- 
ical simulations of a ACDM cosmology, find a dependence of C200 
on both D+ and v = 8 c /a(M, z). They evaluate the overall depen- 
dence of C200 on the linear growth factor, both for their complete 
and relaxed samples, to be D® 5 . The dependence on the linear 
density contrast is also considered to be a power-law, v a . They find 
different values of a for the complete and relaxed samples, namely 
a = —0.35 for the former and a = —0.41 for the latter. Thus, the 
higher the value of D+ the higher the value of the concentration, 
and the lower the value of S c the higher the value of the concen- 
tration. Our findings on the EQp and EQn models do not fit in this 
general picture because both D+ and S c are lower in these models 
than in ACDM, and the differences in S c are less than 1%. 
We hint that the behaviour of EQp and EQn is linked to the redshift 
evolution of the effective gravitational interaction G, as pointed out 
in Sect.f2] In fact, in contrast with ACDM, RP and SUGRA, in EQ 
models the gravitational constant G is substituted by G, which is 
higher (lower) than G at high redshift for positive (negative) values 
of the coupling constant £, while it is equal to G at z — in order 
to recover General Relativity (see Fig. [TJ. This means that in EQp 
gravity becomes weaker at low redshift compared to high redshift, 
while in EQn it becomes stronger. Thus one can expect that in EQp 
(EQn) the halos expand (contract) due to the change in the gravi- 
tational interaction, resulting in lower (higher) concentrations with 
respect to the case in which gravity is constant. 
For the relaxed sample, compared to ACDM, the differences in the 
normalization are less than 10%, while the differences in the slope 
can almost reach 15%, but they are compatible with the standard 
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Table 8. Best-fit parameters, standard errors and reduced chi-squared x 2 of 
the c — M relation equation IIP) for dark matter only density profile fit in 
the region [0.1 — l]/?200 for the complete and relaxed samples of the five 
different cosmological models at z = 0. 



Model 


C8 


A 




B 


0\B 


x 2 


dm (all) 


ACDM 


0.776 


3.59 


0.05 


-0.099 


0.011 


0.48 


RP 


0.746 


3.54 


0.05 


-0.103 


0.011 


1.14 


SUGRA 


0.686 


3.41 


0.05 


-0.098 


0.013 


1.50 


EQp 


0.748 


3.36 


0.05 


-0.097 


0.012 


0.35 


EQn 


0.726 


3.70 


0.05 


-0.069 


0.013 


0.78 


dm (relaxed) 


ACDM 


0.776 


4.09 


0.05 


-0.092 


0.011 


0.66 


RP 


0.746 


4.08 


0.05 


-0.081 


0.011 


0.92 


SUGRA 


0.686 


3.94 


0.06 


-0.081 


0.012 


1.55 


EQp 


0.748 


3.84 


0.05 


-0.097 


0.011 


1.32 


EQn 


0.726 


4.25 


0.06 


-0.081 


0.013 


0.51 



Table 9. Best-fit parameters and standard errors of the c — M relation equa- 
tion (To) for dark matter only density profile fit in the region [0.1 — 1J-R200 
for objects with Af200m > 10 14 Mq h _1 in the complete and relaxed 
samples of the five different cosmological models at z = 0. 



Model 




A 




B 


OB 


dm (all) 


ACDM 


0.776 


3.55 


0.09 


-0.087 


0.038 


RP 


0.746 


3.54 


0.11 


-0.080 


0.046 


SUGRA 


0.686 


3.18 


0.12 


+0.035 


0.060 


EQp 


0.748 


3.35 


0.10 


-0.080 


0.044 


EQn 


0.726 


3.74 


0.11 


-0.114 


0.048 


dm (relaxed) 


ACDM 


0.776 


3.99 


0.10 


-0.055 


0.042 


RP 


0.746 


4.00 


0.12 


-0.017 


0.045 


SUGRA 


0.686 


3.73 


0.13 


+0.051 


0.059 


EQp 


0.748 


3.76 


0.11 


-0.027 


0.046 


EQn 


0.726 


4.19 


0.13 


-0.079 


0.053 
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Figure 13. Ratio between the value of agD + for the ACDM (black), RP 
(blue), EQp (cyan), EQn (red), and SUGRA (green) cosmologies and the 
corresponding value for ACDM as a function of redshift. 

errors. Also in this case, the most extreme cosmologies are EQp 
and EQn, whose normalizations go in the opposite direction with 
respect to their <j$D+, This fact confirms the conclusions we have 
drawn from the complete sample. The values of the reduced chi- 
squared indicate that equation d 1 Ot is a good parametrization of the 
c — M relation for almost all cosmological models. Only SUGRA 
has high values both for the complete and relaxed samples. 
Our results are in good qualitative agreement with the findings of 
iLi et al. ]]), where halos in extended quintessence models have 
lower (higher) concentrations with respect to the ACDM case for 
positive (negative) values of the scalar field coupling. 

We plot the best-fits c — M relations, along with the binned 
data for all the cosmological models, in Fig. [14] We clearly see that 
the results on the normalizations are due to differences in the con- 
centrations over a wide mass range. If we look, for example, at the 
complete sample (left panel of Fig. I14t . we see that the different 
slope of EQn is mainly originated by the less massive bin. But with 
the exception of this bin, EQn shows the highest concentration in 
almost all the mass bins, while in general EQp has the lowest con- 



centration. For the relaxed sample (right panel of Fig. 1141 . the rela- 
tive behaviour of the different cosmological models is even clearer, 
and indeed the differences in the slopes are less pronounced. 

We also try to limit the fit of the c — M relation to the halos 
with K'hoo-m > 10 14 Mq hT 1 , in order to check the effect of in- 
cluding low-mass objects. We report the results in Table [9] For the 
complete sample, we note that the slopes and the normalizations 
are compatible to the fit including also low-mass objects, but that 
the standard errors are a factor of two (four) higher for the normal- 
ization (slope), meaning that the relation is less tight. Notable ex- 
ceptions are EQn, which shows a steep slope, and SUGRA, which 
shows a positive trend of C2oodm with mass (even if with large un- 
certainties) and a consequently low normalization. The slope of 
EQn can be explained by the fact that excluding from the fit the 
lowest mass bin, which is the one with the lowest concentration, 
results in a steeper slope. The behaviour of SUGRA can be ex- 
plained by the fact that this model lacks very massive objects (see 
Paper I), that are the ones with the lower concentration. 

For the same reason, the c — M relation for relaxed objects 
is flatter than when we include low-mass objects. This is expected, 
because we do not have the low-mass objects that have high con- 
centration and we do not have the high-mass objects, which are in 
general more disturbed, that have low concentration. Thus, in this 
mass range, the relation is in general almost flat, with big uncer- 
tainties on the slope. 



7 DARK ENERGY MODELS: RESULTS ON THE TOTAL 
PROFILES 

In this section, we study the impact of baryon physics on the c— M 
relation by analysing the hydrodynamical runs of our simulations 
for all the cosmological models under investigation. This allows us 
to understand the effects of the presence of a dynamical dark energy 
component on the internal matter distribution, including baryons, 
of the halos. We start this analysis with the objects at z = 0, then 
we will consider the redshift evolution of the c — M relation by in- 
cluding also objects at z = 0.5 and z = 1. As we already explained 
in Sect. [3] first of all we fit equation dlO) for the objects at 2 = 0, 
then, keeping fixed the best-fit values of A and B, we fit equation 
dl7t and evaluate the redshift evolution C. 
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Figure 14. (Left panel) The values of C2oodm f° r the complete sample of the ACDM (black), RP (blue), SUGRA (green), EQp (cyan), and EQn (red) 
cosmologies at z = 0. The lines of the corresponding colours are our best-fit of c — M relation equation 1101 and the vertical black bar is the error on the 
normalization of ACDM as listed in Table[7j The symbols in the low part of the panel are the ratios between C2oodm f° r tne model and C2oodm f° r ACDM. 
(Right panel) The same as in the left panel, but for the relaxed sample. 



In TablellOlwe summarize the best-fit parameters, the standard 
errors and the reduced chi-squared of the c — M relations equa- 
tion JlOb and equation J17t for the five cosmological models here 
considered, both for the complete and relaxed samples. For all the 
cosmologies, the values of A are larger than in the dark matter only 
case, indicating that the inclusion of baryons leads to an increase of 
the concentration, while the standard errors remain the same. The 
slope is somewhat flatter than in the dark matter only case for all 
the cosmological models. We already noted both these features in 
the ACDM case (see Sect. [3}: the shallower slope can be explained 
by the fact that the increase in the concentration due to the presence 
of baryons is greater in high-mass objects than in low-mass ones. 
For the complete sample, the total c — M relation reflects the one 
for dark matter, with the normalizations in the same order, with the 
exception of an exchange between EQp and SUGRA. The relaxed 
sample shows higher normalizations than the complete sample, as 
in the dark matter only case, while the slopes are very similar to the 
ones of the complete sample. Also in this case, in comparison with 
the order of <rgD+, the extreme case are EQp and EQn. Thus, the 
trend we find in the dark matter only runs still holds in the hydro- 
dynamical runs. The values of the reduced chi-squared indicate that 
for some models equation dlOt is a good parametrization also for 
the c — M relation of objects extracted from the hydrodynamical 
runs, at least when the complete sample is considered. For relaxed 
objects only SUGRA seems to be described quite well by this rela- 
tion. 

We plot the best-fit c — M relations, along with the binned 
data for all the cosmological models, in Fig. [15] Also in this case, 
we see that differences in the concentration are present over a 
large mass range. For the complete sample, the lower normaliza- 
tion of SUGRA compared to EQp is mainly due to low-mass ob- 
jects, which flatten the relation. The same flattening happens for 
the relaxed sample, leading in this case to a higher normalization 
for SUGRA with respect to RP. 

As for the dark matter only case, we also limit our analysis 
at z = also to objects with M 20 om > 10 14 M hT 1 . We re- 
port the results in Table Qj] We find that the normalizations are 
similar to the one found including also low-mass objects, but the 



Table 11. Best-fit parameters and standard errors of the c — M relation 
equation 4 lOt for total density profile fit in the region [0.1 — lj-Rzou for ob- 
jects with M20OW1 > 10 14 Mq h _1 in the complete and relaxed samples 
of the five different cosmological models at z = 0. 



Model 




A 




B 


0\B 


total (all) 


ACDM 


0.776 


3.74 


0.09 


-0.031 


0.034 


RP 


0.746 


3.67 


0.11 


-0.027 


0.040 


SUGRA 


0.686 


3.49 


0.13 


+0.027 


0.061 


EQp 


0.748 


3.83 


0.10 


-0.130 


0.041 


EQn 


0.726 


3.94 


0.12 


-0.083 


0.045 


total (relaxed) 


ACDM 


0.776 


4.14 


0.10 


-0.007 


0.036 


RP 


0.746 


4.17 


0.11 


-0.004 


0.041 


SUGRA 


0.686 


4.00 


0.15 


+0.069 


0.061 


EQp 


0.748 


4.14 


0.11 


-0.054 


0.043 


EQn 


0.726 


4.41 


0.14 


-0.068 


0.052 



standard errors are higher, and the slopes are generally flatter, with 
some exceptions. The most evident is SUGRA, for which the trend 
of C200 with mass is increasing both in the complete and relaxed 
samples. For the complete sample, EQp and EQn show a steep- 
ening in the slope, a behaviour which is maintained by EQn also 
when considering the relaxed sample. For SUGRA, the explanation 
is the same as in the dark matter only case, i.e. the lack of very high 
mass, low-concentration systems. For EQp and EQn, they show ob- 
jects with low concentration in particular in the low-mass tail, so 
when excluding these objects, the result is a steepening of the slope 
(see Fig.|15t. In all cases, the scatter in the relation is considerably 
higher than when we consider also low-mass objects, up to a factor 
of three (four) in normalization (slope). 

For the hydrodynamical simulations, we also study the red- 
shift evolution of the c — M relation by fitting equation l |17t keep- 
ing the values of A and B fixed to the best-fit values at z = 0. We 
report the results in Table [TO] For all the models, the evolution with 
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Table 10. Best-fit parameters, standard errors and reduced chi-squared \ 2 °f me c ~~ M relation equation (JTOj for total density profile fit in the region 
[0.1 — l]-f?200 f° r me complete and relaxed samples of the five different cosmological models at z = 0, and of the c — M relation equation 117) including 
also objects at z = 0.5 and 2 = 1, keeping the best-fit values of A and B found at z = fixed. 



Model 


0"8 


A 




B 




x 2 


C 




x 2 


total (all) 


ACDM 


0.776 


3.81 


0.05 


-0.061 


0.011 


1.69 


-0.26 


0.02 


2.40 


RP 


0.746 


3.72 


0.05 


-0.073 


0.012 


1.06 


-0.15 


0.02 


1.98 


SUGRA 


0.686 


3.68 


0.06 


-0.057 


0.013 


0.70 


-0.05 


0.02 


3.68 


EQp 


0.748 


3.69 


0.05 


-0.085 


0.012 


0.88 


-0.20 


0.02 


2.99 


EQn 


0.726 


3.94 


0.05 


-0.052 


0.012 


1.55 


-0.21 


0.02 


1.77 


total (relaxed) 


ACDM 


0.776 


4.29 


0.05 


-0.064 


0.011 


1.54 


-0.26 


0.02 


1.97 


RP 


0.746 


4.24 


0.05 


-0.075 


0.010 


2.86 


-0.15 


0.02 


1.84 


SUGRA 


0.686 


4.25 


0.06 


-0.045 


0.011 


1.14 


-0.08 


0.02 


1.11 


EQp 


0.748 


4.13 


0.05 


-0.091 


0.010 


2.02 


-0.16 


0.02 


2.79 


EQn 


0.726 


4.48 


0.05 


-0.057 


0.012 


1.62 


-0.16 


0.02 


1.02 



oil z = reloxed z = 



i 1 1 1 1 — i — i — i i i 1 1 1 1 ■ — -i (—■ — ■— i ■ ■ ■ ■ — ■ — ■ — ■ ■ i 




M 200 ( M s) M 200 ( M o) 

Figure 15. (Left panel) The values of C200tot for the complete sample of the ACDM (black), RP (blue), SUGRA (green), EQp (cyan), and EQn (red) 
cosmologies at z = 0. The lines of the corresponding colours are our best-fit of c — M relation equation 1101 and the vertical black bar is the error on the 
normalization of ACDM as listed in Table|7] The symbols in the low part of the panel are the ratios between C200tot for the model and C200tot for ACDM. 
(Right panel) The same as in the left panel, but for the relaxed sample. 



redshift is negative, meaning that objects at higher redshifts have 
also lower concentrations. For the complete sample, the redshift 
dependence is stronger for ACDM than for the other cosmologi- 
cal models. The two EQ models have similar dependences, lower 
than ACDM, followed by RP and SUGRA, which show a very 
weak evolution with redshift. If we focus on the relaxed sample, 
we see that the ACDM model still shows the strongest redshift de- 
pendence, analogous to the one for the complete sample. RP, EQp 
and EQn have similar evolutions, but while for RP it is the same as 
for the complete sample, for EQ models it is weaker. SUGRA has 
a flatter dependence compared to the other models, but more pro- 
nounced than for the complete sample. The values of the reduced 
chi-squared are rather high for all the cosmological models, both 
for the complete and relaxed samples. Two notable exceptions are 
the values for the relaxed samples of SUGRA and EQn. The results 
shown in TablellOlseem to indicate that equation ( |10t , and its gen- 
eralization equation l !17t . are only a limited parametrization of the 
c — M relation of galaxy clusters extracted from hydrodynamical 



runs of cosmological simulations including dynamical dark energy, 
and of its redshift evolution. 

Finally, we check the evolution with redshift of the normaliza- 
tion A both for the complete and relaxed samples for all the cosmo- 
logical models. We do that by fitting the c — M relation equation 
d 1 0b by keeping fixed the slope at the best-fit value for the com- 
plete sample of ACDM at z = (i.e. B = -0.061, see Table 
[TO) and leaving only A as a free parameter. We compare the dif- 
ferent normalizations found by fitting equation ( llOt in this way for 
both the complete and relaxed samples at z = and at z = 1. 
Thus we can have a snapshot of the imprint on dark energy on the 
concentration of the halos. We summarize the results in Table [T2l 
and plot them in Fig. [16] Indeed we see that, for a given sample 
at z = 0, the normalization is decreasing going from ACDM to 
RP to SUGRA, as expected from agD+. Then, as we already dis- 
cussed, the normalization of EQp is rather suppressed with respect 
to this simple expectation, while the one of EQn is enhanced, due 
to the evolution of the effective gravitational interaction G in these 
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Figure 16. Best-fit normalization comparison for equation d 1 0j with B 
fixed at the best-fit value for the complete sample of ACDM at z = 
for the ACDM (black), RP (blue), SUGRA (green), EQp (cyan), and EQn 
(red) cosmologies. Crosses: total profile fit z = 0, complete sample. Stars: 
total profile fit z = 0, relaxed sample. Triangles: total profile fit z = 1, 
complete sample. Squares: total profile fit z = 1, relaxed sample. The ver- 
tical black bar is the error on the normalization of the complete sample of 
ACDM at z = 0. 



Figure 17. Best-fit normalization comparison for equation d 1 0| with B 
fixed at the best-fit value for the dark matter only complete sample of 
ACDM at z = for the ACDM (black), RP (blue), SUGRA (green), EQp 
(cyan), and EQn (red) cosmologies. Triangles: dark matter profile fit, com- 
plete sample. Squares: dark matter profile fit, relaxed sample. Crosses: total 
profile fit, complete sample. Stars: total profile fit, relaxed sample. The ver- 
tical black bar is the error on the normalization of the dark matter only 
complete sample of ACDM. 



models. If we move instead to z — 1, the relative behaviour of the 
different cosmological models changes. We find that, for both sam- 
ples, the normalization is increasing going from ACDM to RP to 
SUGRA, while EQp is still suppressed and EQn is still enhanced. 
The findings on ACDM, RP and SUGRA at z = 1 are in agree- 
ment with the strong redshift evolution we found for ACDM and 
with the weak redshift evolution we found for SUGRA. There is 
an evolution from a low-normalization to a high-normalization re- 
lation for the first model, and vice versa for the latter. We show the 
values of the reduced chi-squared of the fit as a reference, but we 
do not discuss them because we are imposing the slope for ACDM 
also to other models. 



8 DISCUSSION 

In this section, we discuss the results we have found for the dark 
matter only and hydrodynamical runs for the different cosmologies. 
We start by comparing the normalization of the c — M relation at 
z = for the complete and relaxed samples in both runs. To do this, 
we fix the slope at the best-fit value for the complete sample of the 
dark matter only run for ACDM at z = {i.e. B = —0.099, see Ta- 
ble[8]( and we fit equation < | 1 0I> with only A as a free parameter. We 
report the results in Table [TJl and plot them in Fig. [17] which sum- 
marizes almost all the information on the c — M relation at z — 
for the cosmological models under investigation. Indeed, we see 
that inside each sample, objects in the dark matter only runs have 
a lower normalization than objects in the hydrodynamical runs, in- 
dependently of cosmology and dynamical state. Moreover, inside 
each run, relaxed objects have a higher normalization compared to 
the complete sample. Then, as a general trend, the normalization is 
decreasing going from ACDM to RP to SUGRA, independently of 
the run and the dynamical state. Finally EQn always has the high- 
est normalization while EQp, at least for the dark matter only runs, 
alway has the lowest. In particular, EQn has the highest normal- 
ization also in hydrodynamical runs, while the suppression of the 



concentration in EQp is somehow mitigated in these runs. The be- 
haviour of ACDM, RP and SUGRA is in agreement with the simple 
idea that the normalization of the c — M relation is driven by the 
value of a$D+ (but the one of EQp and EQn is not). We ascribe the 
behaviour of the two EQ models to the different redshift evolution 
of the effective gravitational interaction G, which decreases, going 
from high to low redshift, toward the General Relativity value for 
EQp, while it increases for EQn (see Sect.|2}- This leads to a lower 
concentration, and thus to a lower normalization, for EQn halos and 
to a higher concentration, and thus a higher normalization, for EQp 
halos (see Sect.|6). We show the values of the reduced chi-squared 
of the fit as a reference, but we do not discuss them because we are 
imposing the slope for ACDM also to other models. 

Finally, we compare the concentrations obtained from the dark 
matter only runs with the one obtained from the hydrodynamical 
runs. We do this by taking the ratio between C2ootot and C2oodm in 
each mass bin both for the complete and relaxed samples at z = 0. 
We plot the results in Fig. [18] In the left panel we analyse the ra- 
tio for the complete sample. We see that all cosmological mod- 
els have C2ootot / C2oodm increasing with increasing mass. This fact 
confirms that the baryon physics, as treated in these numerical runs, 
influences more the concentration of high-mass objects compared 
to low-mass ones (see also Fig.[8j. In particular, while ACDM, RP 
and SUGRA have czootot < c.2oadm in some low-mass bins, EQp 
and EQn have in general C2ootot > Caootim, with a less pronounced 
evolution with mass. In the right panel we analyse the ratio for the 
relaxed sample. Here the situation is a bit different. All the cosmo- 
logical models still show a general increase of C2ootot / C2oodm with 
increasing mass, but the evolution is rather shallow. In particular, 
very massive objects in the RP, EQp and EQn models have low val- 
ues of these ratio. In any case, with few exceptions in some mass 
bins, we generally find that C2ootot > 0200dm for all models. This 
analysis indicates that the inclusion of baryon physics in the sim- 
ulations is unable to solve the discrepancy between the predicted 
and observed c — M relation by increasing the slope. Indeed, in 
none of the cases we have analysed, the effect of the batyons is 
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Table 12. Best-fit parameters, standard errors and reduced chi-squared \ 2 of the c — M relation equation (To), with B fixed at the best-fit value for the 
complete sample of ACDM at z = 0, for total density profile fit in the region [0.1 — l]i?200 f° r the complete and relaxed samples of the five different 
cosmological models at z = and 2 = 1. 
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Table 13. Best-fit parameters, standard errors and reduced chi-squared \ 2 of the c — M relation equation (ToJ, with B fixed at the best-fit value for the dark 
matter only complete sample of ACDM at z = 0, for dark matter only and total density profile fit in the region [0.1 — l]-f?200 for the complete and relaxed 
samples of the five different cosmological models at z = 0. 
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to increase the concentration of low-mass objects without affecting 
the one of the high-mass ones. Of course a possible explanation 
of this fact can be that we do not include some kinds of feedback 
in our simulations, in particular AGN feedback. Still, Duf fy et alj 
already showed that none of the different hydrodynamical 
treatments they tried was able to both explain the observed c — M 
relation and the stellar fraction in galaxy clusters. 

Nevertheless, we have seen that different dark energy models 
leave an imprint in the c — M relation. In particular, as already 
noted bv lDolag et all [2004), in ordinary quintessence models, like 
our RP and SUGRA, the normalization of this relation compared to 
the ACDM one is driven by the value of asD+(z co u) with respect 
to the value for ACDM. We find that the same thing is no longer 
valid for extended quintessence models like our EQp and EQn. In- 
deed, in these cases we find that, in presence of a positive cou- 
pling, the value of the normalization is lower than what expected 
from <j$D+ (zcoii) while in presence of a negative coupling it is in- 
creased. We interpret this behaviour as due to the different redshift 
evolution of the effective gravitational interaction G, as explained 
in Sect. [2] Unfortunately, given the discrepancies between the ob- 
served and the predicted c — M relation, it is difficult to use the ob- 
served c— M relation to disentangle different cosmological models 
through the imprints left by dark energy. 



9 SUMMARY 

In this paper we studied the c — M relation for the halos extracted 
from the simulation set introduced in Paper I. From our analysis we 
draw the following conclusions. 

• Concentration: the definition itself of the concentration of a 
halo can lead to very different results. For the reference ACDM 



model at z — 0, we see that the usual fit of a NFW profile, depend- 
ing on the way in which the errors are weighted and the concentra- 
tion is evaluated from the fit, gives values of the concentration with 
differences, at worst, of t he ord er of 5%. The method using V ma x 
discussed in IPrada et al. I J20I 1), instead, returns similar values of 
the concentration only for super-relaxed objects, while in general 
gives higher values of the concentration, up to 20%. The quantity 
J?2500 /R500, which is model independent, is found to be well cor- 
related with the NFW value of the concentration for dark matter 
only simulations, indicating that on average dark matter halos are 
indeed well described by a NFW profile. In a given mass bin, the 
mean concentration is higher and the intrinsic scatter is lower for 
objects that are in a relaxed dynamical state and well described by a 
NFW profile. Halos extracted from the hydrodynamical runs have, 
in general, higher concentrations than their dark matter only coun- 
terparts. The effect is more evident in high-mass systems, due to a 
higher baryon fraction in the inner regions. 

• c — M relation: limiting our analysis to the ACDM model, 
there is a large intrinsic scatter in the values of the concentration 
for objects of a given mass, ranging from 15% up to 40% depend- 
ing on mass and dynamical state. Nevertheless, if we consider the 
mean value of concentration in a given mass range, there is a good 
correlation between concentration and mass. The concentration is 
slightly decreasing with increasing mass, and this relation is well 
described by a power-law, with two free parameters, the normaliza- 
tion and the slope. The normalization, that we define as the value of 
the concentration of a halo with a mass of 10 14 M©, is increasing 
going from A — 3.59 for the complete sample to 4.09 and 4.52 
in the relaxed and super-relaxed samples, respectively. The slope 
is similar for the different samples considered and is B ~ —0.1. 
In general, we find a good agreement with the results from other 
works in literature. The normalization is higher by 5 — 15% for 
objects extracted from the hydrodynamical runs compared to the 
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Figure 18. (Left panel) Ratio between C2ootot and c 2 oodm f° r the complete sample of the ACDM (black), RP (blue), SUGRA (green), EQp (cyan), and EQn 
(red) cosmologies at z = 0. (Right panel) The same as in the left panel, but for the relaxed sample. 



one of the dark matter only runs. The slope is flatter by 30% for 
the hydrodynamical simulations than for the dark matter only ones. 
This is expected because massive objects from the hydrodynamical 
simulations have higher concentrations than their dark matter only 
counterparts. 

• Redshift dependence: the c — M relation shows an evolution 
with redshift, with concentration decreasing with increasing red- 
shift as (1 + z) -0 ' 3 . For the ACDM model, considering also objects 
at 2 = 0.5 and z = 1, the redshift evolution is more pronounced 
for the dark matter only simulations than for the hydrodynamical 
ones, while it is similar for the complete and relaxed samples. 

• Dark energy models: we find that the normalization of the 
c — M relation in dynamical dark energy cosmologies is different 
with respect to the ACDM one, while the slopes are more compat- 
ible. In particular, at z = 0, the differences in the normalizations 
for RP and SUGRA when compared to ACDM reflect the differ- 
ences in agD+, with models having a higher agD+ also having a 
higher normalization. This simple scheme is not valid for the EQp 
and EQn scenarios. In the former case, the normalization is lower 
than expected considering <j$D+, while in the latter it is higher, and 
indeed EQn is always the model with the highest normalization, re- 
gardless of the dynamical state of the objects or the runs they are 
extracted from. This behaviour is due to the different redshift evo- 
lution of the effective gravitational interaction G in these models. 
Indeed, going from high to low redshift, G decreases (increases) for 
EQp (EQn) toward the value for ACDM and ordinary quintessence 
models, thus decreasing (increasing) the concentration of the halos. 
This is a very important result because it shows a direct manifes- 
tation of the coupling between gravity and the quintessence scalar 
field. For objects extracted from the hydrodynamical runs, we also 
study the redshift evolution of the c — M relation. We find different 
evolution for different dark energy models, with values of C rang- 
ing between —0.05 and —0.26. ACDM has the strongest evolution, 
while SUGRA has the weakest, while RP lies in between. For these 
three models, the behaviour is similar for the complete and relaxed 
samples. EQp and EQn models, instead, show an evolution simi- 
lar to ACDM for the complete sample and similar to RP for the 
relaxed sample. It is interesting to note that at z — the normal- 
ization decreases from ACDM to RP to SUGRA, while at z = 1 
the situation is completely reversed, indicating a stronger redshift 
evolution in SUGRA with respect to ACDM. Independently of red- 



shift, EQp has always the lowest normalization while EQn has al- 
ways the highest, meaning that the effect of the coupling between 
gravity and the quintessence scalar field is important also at higher 
redshift. 
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